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The formation of spherical polaron clusters is studied within the Frohlich polaron theory. In a 
dilute polaron gas, using the non-local statistical approach and the polaron pair interaction obtained 
within the Pekar strong coupling theory, the homogeneous phase results to be unstable toward the 
appearance of polaron clusters. The physical conditions of formation for the clusters are determined 
calculating the critical values of electron-phonon interaction for which bound states in the collective 
polaron potential develop. Finally the sequence in the filling of the states is found and the stability 
of the clusters is assessed. 
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I. INTRODUCTION 



In the last years the presence of strong electron-phonon coupUng and polaronic effects has been pointed out by 

many experimental results in several compounds, such as high-tempcrature cuprate superconductors.^ colossal magnc- 
toresistance manganites,^ nickelates and quasi-lD materials organic conjugated polymers.^ Furthermore it has been 
debated if the electron-phonon interaction can give rise in perovskite oxides to charge-ordered states or to more 
complex electronic phases such as stripes, strings or clustered states.^ 

The large amount of experimental data has renewed the interest in studying simplified electron-phonon coupled 
systems of the Holstein^ or Frohlich^ type. The formation of the polaron have been long studied in the frame of the 
Frohlich model, where the polar long-range interaction between electronic and ionic charges is taken into account and 
the medium is considered continuous. For its relation to theories regarding cuprate superconductors, large attention 
has been devoted to the formation of the bipolaron (bound state of two polarons), establishing the range of the values 
of the electron-phonon coupling constant and of the dielectric parameters of medium which allow its existence. Also 
the binding energy, the effective mass, its internal structure and optical features have been studied. Actually 
many methods have been used to study the bipolaron formation in the Frohlich scheme: Lce-Low-Pines approach for 
intermediate values of the electron-phonon coupling constant,^ Pekar polaron strong coupling theory,^'^^'^^ and path- 
integral technique. They give the possibility to define an effective polaron-polaron potential, which is repulsive both 
at large and small distance between the particles, while it is attractive at intermediate distance. This polaron-polaron 
interaction has been used as starting point for many-body calculations finding that, for densities smaller than those 
typical of metals, the many-polaron system exhibits a charge density instability in the intermediate coupling regime. 

A related research activity concerns the study of the conditions under which electronic or polaronic Wigner 
crystals, ^^'^^ polaron molecules or clusters^^ and strings^^ can form. Indeed in the regime of very low densities 
the Wigner crystal of electrons can transform into a polaronic crystal by increasing the electron-phonon coupling. "'^^ 
However increasing the density at strong electron-phonon coupling the Wigner crystal of polarons becomes unstable. 
Such an instability suggests that novel types of electronic or polaronic structures such as clusters, molecules or strings 
can be stabilized by a strong electron-phonon interaction. Actually the molecules can arise at conditions determined 
by the structure of an effective interaction at short and long distances.^® Their behavior is very quantum, and they 
can assume a stringed form. It has been suggested that in perovskite oxides, such as manganites and cuprates, these 
molecules could represent the intermediate step of more complex structures such as stripes. "'^^ Finally the electronic 
and ionic structure of Metal-Ammonia solutions^*'^^ represents a system where polaronic clusters can be actually 
observed. 

In this work we start from the knowledge of the effective polaron-polaron interaction within the Pekar strong 
coupling theory.^'^'^^'^^ Using the non-local statistical approach, ^^'^^ we show that in the homogeneous phase the 
mode due to the collective polaron potential is unstable suggesting that polaron clusters can form in the system. 
Next we show that there is a range of values of the electron-phonon interaction and of the static and high frequency 
dielectric constants which allow the formation of clusters with a number of polarons larger than two. Finally we 
discuss some features of the cluster such as its radius and its stability. 

In sect. II we discuss the basic equations; in sect. Ill we analyze the collective excitations of the homogeneous 
polaron phase and its instability due to the attractive interaction between polarons; in sect. IV we indicate the 
conditions for which the polaron cluster forms, we discuss the shell structure in the cluster and finally we study the 
stability of the polaron clusters. 



II. BASIC EQUATIONS 

In the system we have two types of particles, the polarons and the ions, with distribution functions fp{f,v,t) and 
fi{f,v,t), respectively, in the single-particle phase space. In the non-local and non-linear Vlasov kinetic equations^^ 
the distribution functions fp{f, v, t) and fi{f, v, t) are solutions of the equations 

-^^^ ' ' ^ + V ■ Vrfp{f, V, t) V. [[/(r, t) - e$(r, t)\ ■ V./p(r, v,t) = Q (1) 



dMr v,t) v./,(f,i;,i) - — V.$(f,t) ■ V./,(f,tT,t) = 0, (2) 

Ot rrii 

where rup and rrii are the effective masses of polaron and ion, respectively, e is the electron charge, and U{f, t) is the 
collective self-consistent potential 
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U{f,t) = j dr'dv'K{\f- 



■r-'|)/p(r^tf,t), 



(3) 



with the kernel K indicating the polaron-polaron potential depending only on the relative distance between particles. 
The distribution function fp is determined not only by the polaron collective potential U but also by the electrostatic 
potential $ arising from the non- punctual compensation between the polaron and ionic charges. Indeed the potential 
$ satisfies the Poisson equation 



Aire , ^ , 




where the density p(f, €) is defined as 



p(f, t) = pi{r, t) - pp{r, t)= J dvfi{r, '"^t) - j 



dvfp{f,v,t). 



(4) 



(5) 



We note that the force acting on the ions is related only to the electrostatic potential $ screened by the static dielectric 
constant e^. 

In order to solve the Vlasov equations, the polaron-polaron potential K has to bo specified. We use the potential K 
obtained within the Pekar strong coupling bipolaron thcory®'^^'-'^^'^^ which has been adapted in order to find reliable 
results in the intermediate to strong coupling regime. This polaron-polaron potential is attractive at intermediate 
distances, this behaviour resulting by the Prohlich interaction and effects of quantum mechanical exchange and inter- 
electron correlations. Moreover the potential depends on the ratio e*/eoo, where e* is defined through the equation 
1/e* — I/eoo — I/es, with e^o high frequency dielectric constant. In Fig. 1 we show the potential K for different 
values of the ratio e* /e^o when the bipolaron is in its singlet ground-state and the asymptotic energy of the two free 
polarons has been subtracted. Throughout the paper, the values 1.10, 1.08, 1.05, 1.02 and 1.00 are used for the ratio 
e*/eoo, i.e. the ratio Cs/eoo varies from 10 to infinity. 

The potential K obtained within the Pekar strong coupling theory can be approximated by the following analytical 
expression 



K{r) = 2a^nuJo 



(e*/eoo - l)aS 



-Sr 



where Kx (r) is given by 



Kr{r) = d + cr^ + [a + h{r - pof] (l - e"^'') 



(e*/ecx, - l)al 



(6) 



(7) 



In Eq. (6) a = (l/2e*)(e^/?ia;o)(2TO*a;o/?i)^'''^ is the dimensionless electron-phonon coupling constant, m* is the 
effective mass at the bottom of the conduction band, is the longitudinal optical phonon frequency in the long 
wave-length approximation, and Oq = Q!~^(?i/2m*cc;o)^/^ is the effective Bohr radius. The parameters a, h, c, d, 7, 5, 
and po of the potentials in Eq. (6) and (7) have been obtained by an accurate fit of the polaron-polaron potential in 
the Pekar theory.^'^^'^^'^^and are listed in Table I. In the numerical calculations it is furthermore assumed e^o = 2, 
m* — m with m electron mass, and hujo = O.Q3eV. 

The Vlasov equations in Eqs. (1) and (2) are valid in the classical limit and the use of polaron-polaron interaction 
(6) is accurate for a dilute system. Therefore the temperature T has to be high compared with the degeneracy 
temperature of polarons: T > %^Np^^ /nipkB, with Np polaron density and ks Boltzmann constant. This condition 
is verified at reasonable temperatures if A^p < 5 x lO^^cm"^. We note that for these densities some energy scales, for 
example the polaron plasmon frequency, are smaller or of the same order than the optical phonon frequency. Clearly, 
for the temperatures considered above, quantum effects in the statistics can be neglected and the Boltzmann-Maxwell 
distribution function can be used at the thermodynamic equilibrium. 



III. COLLECTIVE EXCITATIONS OF THE POLARON SYSTEM 



In this section we analyze the collective excitations of the polaron system. As discussed in the following subsection 
A, the polaron plasmon is obtained neglecting the role of the collective potential U in the Vlasov equations (1) and 
(2). On the contrary in subsection B we will examine the collective excitations induced by the potential U in absence 
of the electrostatic potential These density oscillations are completely different from the polaron plasmons, since 
at small values of the momentum k they have a dispersion relation approximately proportional to k. 
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The procedure of calculation of these two types of excitations is the same. Indeed only the role of the polarons is 
mainly considered and the Vlasov equations arc linearized taking 



fp{r,v,t)=fo{v)+^{r,v,t), (8) 
where fo{v) is the Maxwell distribution function 

associated to the homogeneous spatial distribution of the particles and (j){f, v, t) is strongly smaller than fo{v). There- 
fore it is assumed that the perturbations induced by the time-dependent potentials are weak. 



A. Polaron plasmons 

The polaron plasmons are the collective excitations in the limit of collective potential U zero. In this case the external 
force affecting the polaron distribution can be simply related to an electrical field ^(r, t) defined as £^(r, t) = — V$(r, t). 
Neglecting second order terms, the equation (1) becomes 

+tT. VrcP{r,v,t) = —E{f,t) ■ V./o(tT). (10) 
at rup 

Fourier transforming in space and time, one obtains 

<f>ik,v,w) = ^^j^""^ ^ ■ VJoiv), (11) 
imp[k ■ V — Lu) 

where the coefficient of V^,/o(^/) represents the amplitude of momentum that the polaron acquires under the field 
E. Clearly this quantity has to be small compared with the mean momentum obtained through the equilibrium 
distribution fo{v). 

Using Eq. (5), the Fourier transform of particle density p can be evaluated. This last quantity is connected to 
the polarization P of the system since ik ■ P = —j-p- Employing the pole-rule defined by Landau,'^* we obtain the 
relation 

ik-P{k,Lo) = -—E{k,Lo)-[dv YjlM^l , (12) 

J imp{k ■ V — (j) —irj 

where t] is infinitesimal. If the field E and P are directed along k, then the longitudinal dielectric constant ei{k,u}) 
can be derived by means of the equation A-kP = {ei — l)E, yielding 

ei{k,uj) = l--—^ dv- ^ -. (13) 

^si^ J imp{k ■ V — CO) — IT] 

The dielectric constant of Eq. (13) is a complex quantity implying that the energy of the electrical field can be 
dissipated in the medium (Landau damping). Hence longitudinal electrical waves can propagate through the system 

and their dispersion relation is obtained by the equation e/(fc, w) = 0. Supposing u » kvx, with Vt = {ksT /mpY/"^ 
mean quadratic velocity along the direction of propagation, the zeros of the equation can be derived. At low order in 
momentum the real part of the frequency is 

uj = Upp{l + ^k''al), (14) 
with Up = [eskBT/AnNpe^y/^ and cOpp polaron plasma frequency defined as 
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while the imaginary part is exponentially small. As expected in the limit ojpp < ujq, the plasma frequency is screened 
by the static dielectric constant e^.^^ 

In this derivation wc have taken into account only the polaron contribution. Within the same approach it would 
be possible to consider also the role of the ionic plasmons through Eq. (2) and the interplay between polaronic and 
ionic oscillations.^^ However, it is more interesting to focus on the collective excitations due to the polaron collective 
potential U. 



B. Excitations due to the collective potential U 



In this subsection we consider the time-dependent kinetic equations (1) for the polaron distribution function ne- 
glecting the contribution from the electrostatic potential $. We investigate the propagation of longitudinal waves 
with lengths larger than the polaron-polaron distance in the bipolaron. Therefore our aim will be the calculation of 
the dispersion relationship linking frequency w and wave vector k in the limit A; — > 0. 

If one neglects second order terms, the Vlasov equation (1) becomes 



d(t){f, V, t) 
dt 



1 



V„/o(^/)-V, / dr^dv'Ki\r-r^\)cl)ir^,v',t). 
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The solution of the integro- differential equation (16) can be written in the following form' 
Hr,v,t) = -^Vyfoiv)-Vr [ d7^\ [ dTK{\r-v{t-T)-r^\)p{r^,T)+ f dTK{\f - v{t - t) - 7^\)p{f^ ,t) 

^"^P J Uto Jtl 

with to <t <ti and p{f, t) defined in Eq. (5). Integrating with respect to the velocity, one obtains 



(16) 



, (17) 



P{r,t) 



^ [ diNyfoiv) -Vr [ df [ dTK{\r-v{t-T)-r'\)p{f',T)+ f dTK{\r - v{t - t) - r'\)p{r' ,t) 

^"^P J J iJto Jtl 



Considering the spatial and temporal Fourier transform of the density p, Eq. (18) becomes 

p{k, w) = - / rfTe-''^(*-^)G(A;, t - T)p{k, cj) + - dTe-''^'-*-''^G{k, t - T)p{k, w), 
where G{k, t) is defined by 

G(fc, t) 



^fc. f dve'^-^'VMv), 
mp J 



(18) 



(19) 



(20) 



with a{k) Fourier transform of the polaron-polaron potential shown in Fig. 2 for different values of the ratio e*/eo 
The condition of existence of the nontrivial solutions of the equation (19) is the following 



POO 

/ dtG{k,t)cos{ujt) = 1, 
Jo 



that links implicitly the frequency lo to the momentum k. We obtain the dispersion for small k in the form 



/ V' 

" V2T 



vt 

76 



k 1 



keT 

NpGik) 



1/3 



(21) 



(22) 



where vt = {ksT /mp)^/'^ is the mean quadratic velocity along the direction of propagation. Therefore the system 
exhibits a collective mode sustained by the coherent self-consistent interaction arising from neighboring particles. 
Clearly these excitations are well defined if 



ksT 



(23) 
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As deduced from Fig. 2, the collective mode can certainly propagate for values of k around where cr(fc) is positive. 
The Fourier transform a at fc = is given by the spatial integral of the potential K that, as shown in Fig. 1, 
becomes completely negative as the ratio e*/eoo approaches the unity. By fixing the ratio e*/eoo and increasing the 
electron-phonon coupling constant a, the potential a{k) deepens implying that the region of stability in k is reduced. 
So there is a range of values of the momentum where Eq. (23) is violated. This suggests the possibility that the 
particles in equilibrium are not longer single polarons, but, due to the attractive effect of the interaction, a polaron of 
the system can bind another one in order to form a bipolaron. Clearly two particles could form with a third a cluster 
of more than two particles. Therefore it is important to study directly the formation of polaron clusters induced by 
the collective potential U. 



IV. POLARON CLUSTERS 



In this section we will study the solutions of the Vlasov equations (1) and (2) in the temporal stationary regime. 

Clearly these equations are satisfied if the polaron fp and ion fi distribution functions are independent of position. 
Our aim is to study the solutions of equations (1) and (2) near to the homogeneous ones. 
In stationary conditions, we write for the polaron and ion distribution function 



(24) 



where j stands for p or i, pj{r) and Wj{v) the spatial and velocity distribution functions, respectively. If the velocity 
distribution function is that of the equilibrium, then from equations (1) and (2) the spatial distribution functions of 
polarons and ions can be deduced as 



and 



Pi{r) = Nit 



(25) 



(26) 



where Ni is the ion concentration, such that Np = Ni = Nq since the system has charge neutrality. Substituting (25) 



and (26) in (3) and (4), we obtain the coupled equations 

U{r) = j df'K{\f-f'\)pp{f') 



and 



A$(rO = 



_ / 47re\ 



[Pp{r) - Pi{r)] ^ 



(27) 



(28) 



that allow to determine self-consistently both the collective polaron U{f) and electrostatic "I>(r) potentials. 

The cixact solution of Eqs. (27) and (28) is a formidable task. Therefore, considering the dependence of the potential 
K on the relative distance, we iteratively build a solution for a spherical symmetric cluster of radius R obtaining at 
the lowest order ^o{r) = and 



U, 



°^^'~\NoJdr'K{\f^-r\)\, r > R. 



(29) 



We note that the collective potential Uo{r) is directly proportional to the particle density A^o and through K it 
depends on the electron-phonon coupling constant a. Through the use of the coordinate transformation 



+ — 2xr cos(6'), sds = xr sm[6)d6, 



the potential inside the sphere C/q™^ (r) is written as 



^0 (»") = — — 



pr rr+x 

/ dxx I ( 

Jo Jr-x 



dssK{s) 



dxx 



dssK{s) 



(30) 



(31) 
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and connects with continuity to that outside [7'^°"*^ (r) 

U^°^'\r) = r dxx ['^\ssK{s),r > R. (32) 

T Jo Jr-x 

The local modifications of the spatial distribution function due to the attractive forces are described by the polaron 
spatial density A'^i (r) 

= (-if) (^') 

showing that the distribution of the polarons is determined by the collective potential UQ{r). Actually the charge 
neutrality is locally perturbed and the increase of the polaron density can be consistent with a bound cluster of 
particles. 

In Figs. 3(a) and 3(b) we show the collective potential for different radii R and for the value of e*/eoo = 1-05. The 
qualitative behaviours are similar in all the considered cases except for e* /coo = 1-00. For a fixed ratio €* /^oo, Ua{r) 
is negative in a region around r = 0, it assumes the minimum value at r = and becomes positive at large values of r. 
The quantity C/o(0) is first decreasing and then increasing as function of R, so that it exists a value of R which gives 
the minimum value of Uq{Q). Only in the case e*/eoo = 1 (i-e. Cg = oo) Uo{0) is always decreasing as function of R. 
In Fig. 4 we show the quantity Uq{0) as function of R for all the indicated values of e* /e^o, except for e*/eoo = 1- It 
occurs that the decrease of t* j too implies the increase of R. We give to R the meaning of the best dimension of the 
cluster. 

It is important to calculate the minimum values of a for which the cluster begins to form. We have evaluated 
the critical values ac within the full quantum mechanical approach of a particle in the collective potential UQ(r) 
calculating when the first bound level develops.'^" Hence the quantity ac indicates the minimum value for which the 
bipolaron forms. It has been checked that these values are only slightly smaller than those obtained in the next 
subsection through the statistical model that can be generalized at arbitrary levels. In particular we will investigate 
the resulting shell structure in the cluster. 



A. Shell structure of the poleiron cluster 



In analogy with the structure of nuclei^^ and metal clusters, we try now to calculate the conditions for existence 

of the polaron clusters and the related shell structure. For the spherically symmetric collective potential ljQ(r\ we 
use the same criterion of the atomic physics regarding the quantum numbers to introduce in order to characterize the 
state of a single electron in the atoms. In other words we introduce the quantum numbers of single particle states 
without considering neither the electrostatic corrections nor the relativistic ones. In the statistical model,^^'^'' the 
shells with the orbital quantum number i (/ = 0, 1, 2, 3...) begin to form when Ni defined by 



Ni = 



4 

7r?i 









/ dr 




J Ri 



-2mpUo{r) - 













[ dr 






J Ri 



(34) 



becomes larger than 1. In the above formula R is the cluster radius, i?i is the lower value of r for which the function to 
integrate is zero, P is the maximum polaron momentum. In the statistical approximation the maximum momentum 
P is connected with the local polaron density (r) given in (33) through the relation 



P(r) = 27r?i 



3iVi(r) 
Stt 



1/3 



(35) 



so that, writing Ni{r) = mpo{r), where m is the number of the particles involved in the cluster, we have to ensure 
that J dfpo{r) = 1. Taking into account Eq. (35), we find the following equation 



1 Id 
r 3dr 



(36) 



which allows to determine explicitly the value 

In table II we have reported the calculated values of the electron - phonon coupling constants a^; from which the 
shell with the indicated value of I begins to be filled. We find that, increasing the polaron concentration, the collective 
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potential Uo{r) deepens, so that the values of ad decrease. For example, taking e*/eoo = 1.05, for Nq = 10 cm^'' the 
first shell I = starts to be filled as soon as q;c,o = 10.2, while for Nq = lO^^cm"*^ we have a^o = 18.2. Furthermore, 
since the increase of the number I reflects a higher binding energy of the particles in the cluster, in table II the values 
of ac,i get enhanced as a function of I. 

Finally we can calculate the sequence of the shell formation in a polaron cluster. In analogy to the filling of states 
for nuclear matter and metal clusters, ^^'^^ we obtain the sequence (ls)^(2p)^(3(i)^°(2s)^(4/)^'* for the polaron cluster. 
This order differs from the known order of filling of electron levels in atoms: in fact shells with 2, 8, 18, 20 polarons 
are stable. Then it is possible that in the system stable clusters can form becoming the basic units in equilibrium at 
a given temperature. 

Finally we can make a rough estimate of T^ri the highest temperature under which the clusters can form or, 
better, the lowest temperature above which the polaron phase is homogeneous. Actually, knowing now the critical 
electron-phonon coupling constants that trigger the cluster formation, we can calculate the temperatures for which 
the collective mode discussed in the previous section becomes unstable signaling the possibility of polaron clusters. 
The stability condition of the collective mode of the homogeneous phase is violated when a{k) vanishes becoming 
negative at larger values of k. Therefore the polarons are the carriers at equilibrium approximately for temperatures 

T>T,. = *&^, (37) 

Kg 

where \(j{k)\min is the modulus of the Fourier transform of the polaron-polaron potential K calculated for the value 
of k corresponding to its minimum value (the attractive part drives the instability) and for the critical values ac,i of 
the coupling constant. In the table III the critical temperatures Tcr are shown. We note that the critical temperature 
obtained in Eq. (37) is a function of the polaron density also through the constants ac,i that, as reported in tab. II, 
are strongly dependent on A^o- Since the collective potential responsible for the cluster formation is proportional to 
A'o, we expect that the critical temperatures strongly increase as a function of the polaron concentration, as confirmed 
by the results of tab. III. Clearly it is possible to define the critical temperatures for condensation of polarons for 
different shells. Increasing the quantum number I, the binding energies becomes larger, the cluster becomes more 
stable, so that the critical temperatures are enhanced. 

In this section we have realized that, at the lowest order of the iterative procedure used for solving equations (27) 
and (28), the electrostatic potential $ defined in Eq. (4) does not affect the cluster formation. However, at higher 
orders, the difference between the density of polarons and ions could perturb the states of the cluster. The frequency 
scale of such perturbation is close to the frequency Wpp of the polaron plasma oscillation given in Eq. (15). Such 
fluctuations can contribute to the broadening of the states determined in the cluster. Obviously the broadening is 
small if ijjpp < LOi, where coi is the frequency of motion of the particle in the collective potential at the level with a 
quantum number I. From this inequality it is possible to find the highest value of the electron - phonon coupling 
constant a, such that for a < a the cluster is unstable against these perturbations. For example, at the concentrations 
A^o = 10^"^ cm~^ and Nq = lO^^cm"^ for / = and e*/eoo = 1,08, we have a = 16.6 and a = 12.5, respectively. These 
results do not contradict the data listed in the table II. Therefore when the polaron cluster is formed in the system, 
it is also stable against the energy fluctuations of the bound states of the collective potential Uq- 



V. CONCLUSION 



Using the statistical model, we have found that clusters can form in a dilute polaron gas. We have calculated the 
critical values of the electron-phonon interaction for which the ground and excited state shells can be filled. The found 
sequence of the shells is analogous to that of the nuclei and of the metal clusters rather than that of the electron 
states in atoms. We have also calculated the stability of the polaron clusters against the fluctuations of the collective 
potential. We have found that the values for which the cluster can form are sufficiently large to assure the stability of 
the system. Finally the critical temperature for the condensation of the polarons in the cluster has been calculated. 
Our calculation scheme based on the strong coupling Frohlich electron-phonon theory is coherent with the results 
since all critical values for the formation of the clusters are found in the regime of large coupling. Clearly this implies 
that the theoretical results can be discussed in connection with real systems only in situations of poorly screened 
strong electron-phonon coupling. ^^■^'^ 

The system with polaron clusters is characterized by collective modes different from those of the homogeneous phase, 
for example that induced by perturbations which tend to deform the cluster. In such a case the difficulty is to solve 
the time dependent equations starting at the zero order from the spatially perturbed polaron distribution function 
calculated in the previous section. Actually vibronic modes have been discussed in electronic bound configurations 
showing that their frequency is of the order 10^^ — lO^^Hz for typical parameters of metallic oxides.^® Finally we note 
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that, in the quantum case of high density and low temperature characteristic of perovskite oxides, the system can 
exhibit a true phase transition toward a clustered state that in the one-dimensional case can be a string. 



FIGURE CAPTIONS 

Fig.l Polaron-polaron potential K (in units of 2a^h,ujQ) as a function of the distance p (in units of Oq) for different 
values of the ratio e* /eoo- solid {e* /too = 1-00), large dot {e* /too = 1-02), dash (e*/eoo = 1.05), dash-dot 
(e7e<x> = 1.08), small dot (e7e<x> = 1.10). 

Fig. 2 Fourier transform a (in units of 20^^000 Og^) of the polaron-polaron potential if as a function of the momentum 
k (in units of flg"^) for different values of the ratio e* jtoo'- 1-soHd {e* /e^o = 1.10), 2-dot {e* /too = 1.08), 3-dash- 
dot (e*/eoo = 1-05), 4-dash = 1.02), 5-solid = 1.00). 

Fig. 3 (a) The collective polaron potential Ua (in units of 2a^hLUQ/Na) as a function of the distance r (in units of Aq) 
for e*/eoo = 1.05: 1-solid (radius of the cluster R = 7«q), 2-dash (radius of the cluster R = Hag), 3-dash-dot 
(radius of the cluster R = 15ag). Optimal radius of the cluster is 15ag. The dot curve defines the boundaries of 
the cluster. 

(b) The collective polaron potential Uq (in units of 2a^fiu!o/NQ) as a function of the distance r (in units of ag) 

for e* jf-oc = 1.05: 1-solid (radius of the cluster R = 16ag), 2-dash (radius of the cluster R = 20ag), 3-dash-dot 
(radius of the cluster R = 22ag). Optimal radius of the cluster is ISajJ. The dot curve defines the boundaries of 
the cluster. 

Fig. 4 The collective polaron potential Uq (in units of 2a^huJo/NQ) at 7' = as a function of the cluster radius R (in 
units of Uq) for different values of the ratio e* /coo- 1-dot (e*/eoo = 1-10), 2-dash-dot (e*/eoo = 1.08), 3-dash 
(eVeoo = 1.05), 4-solid {e*/eoo = 1.02). 
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TABLE I. Parameters of the polaron-polaron potential K 





a 


6 (ar^) 


c (a* ^) 


d 


7 


5 {a*r') 


PQ 


1.10 


-7.95 X 10"' 


6.875 X 10"^ 


-6.821 X 10"'' 


1.4 X 10"-' 


0.2475 


0.2047 


7.95 X lO"'' 


1.08 


-5.175 X lO"' 


6.885 X 10"'' 


-6.857 X 10"^ 


0.001 


0.245 


0.20 


8.17 X 10"^ 


1.05 


-1.907 X 10-" 


6.855 X lO-' 


-6.852 X 10-" 


8.583 X 10"* 


0.25 


0.19 


8.81 X lO"'' 


1.02 


-6.0 X 10"" 


6.9188 X lO-" 


-6.887 X 10-" 


-3.097 X lO-' 


0.286 


0.166 


9.25 X lO"'' 


1.00 


-7.92 X Kr* 


6.709 X lO"-' 


-().()77 X lO"'' 


-7.750 X 10"- 


0.:!2 


O.l.") 


0.1087 



TABLE IL The critical values of the electron-phonon coupling constant when the shell with momentum I is filled 



t / < JO 


Ek'ctroii-plioiiou (XMipliiig const aiit tt, .; 


No = lO^^cm"^ 


I = 


I = 1 


/ = 2 


I = 3 


1.10 


76.0 


89.8 


110.7 


129.4 


1.08 


51.5 


59.8 


73.7 


86.1 


1.05 


32.2 


35.0 


42.2 


48.9 


1.02 


23.7 


25.4 


30.4 


35.1 


1.00 


17.5 


18.3 


21.6 


24.7 


A^o = W"cm-' 




1.10 


42.7 


50.5 


62.3 


72.7 


1.08 


29.0 


33.6 


41.5 


48.4 


1.05 


18.2 


19.7 


23.7 


27.5 


1.02 


13.3 


14.3 


17.1 


19.7 


1.00 


9.8 


10.3 


12.2 


13.9 


Ao = 10'** cm-^ 




1.10 


24.0 


28.4 


35.0 


40.9 


1.08 


16.3 


18.9 


23.3 


27.2 


1.05 


10.2 


11.1 


13.4 


15.5 


1.02 


7.5 


8.0 


9.6 


11.1 


1.00 


5.6 


5.8 


6.8 


7.8 



T^VBLE 111. Critical tcnipt-rattircs (in niiits of 2n ti^ii/kn) corresponding to tlic filling of tlit^ dnstcr .slidl witli monicntiini /. 





Critical temperature Tc, 


Ao = lO'" cm"^ 


/ = 


/ = 1 


I = 2 


I = 3 


1.10 


2.4 


3.1 


4.7 


6.5 


1.08 


1.6 


2.0 


3.0 


4.1 


1.05 


0.91 


1.1 


1.6 


2.1 


1.02 


0.72 


0.83 


0.6 


0.8 


1.00 


0.84 


0.92 


1.2 


1.7 


Ao = 10" 'cm-" 




1.10 


7.0 


9.8 


15.0 


28.4 


1.08 


4.9 


6.2 


9.5 


12.9 


1.05 


2.9 


3.4 


4.9 


6.7 


1.02 


2.3 


2.6 


3.8 


5.0 


1.00 


2.6 


2.9 


4.1 


5.3 


Ao = lO'** cm"^ 




1.10 


22.2 


31.1 


47.2 


64.4 


1.08 


15.6 


19.6 


29.8 


40.6 


1.05 


9.2 


10.8 


15.8 


21.1 


1.02 


7.2 


8.2 


11.8 


15.8 


1.00 


8.6 


9.3 


12.7 


16.7 
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